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If metal surfaces are exposed to sufficiently oxygen-rich environments, oxides start to form. Impor- 
tant steps in this process are the dissociative adsorption of oxygen at the surface, the incorporation 
of O atoms into the surface, the formation of a thin oxidic overlayer and the growth of the once 
formed oxide film. For the oxidation process of late transition metals (TM), recent experimental 
and theoretical studies provided a quite intriguing atomic-scale insight: The formed initial oxidic 
overlayers are not merely few atomic-layer thin versions of the known bulk oxides, but can exhibit 
structural and electronic properties that are quite distinct to the surfaces of both the correspond- 
ing bulk metals and bulk oxides. If such nanometer or even sub-nanometer surface oxide films are 
stabilized in applications, new functionalities not scalable from the known bulk materials could cor- 
respondingly arise. This can be particularly important for oxidation catalysis, where technologically 
relevant gas phase conditions are typically quite oxygen-rich. In such environments surface oxides 
may even form naturally in the induction period, and actuate then the reactive steady-state behavior 
that has traditionally been ascribed to the metal substrates. Corresponding aspects are reviewed 
by focusing on recent progress in the modeling and understanding of the oxidation behavior of late 
TMs, using particularly the late Ad series from Ru to Ag to discuss trends. 



I. INTRODUCTION 



Confinement to nanoscale dimensions and the result- 
ing quantization of electrons gives rise to two-dimensional 
quantum films, one-dimensional quantum wires, or zero- 
dimensional quantum dots. All three situations may ex- 
hibit properties and functionalities that make them dif- 
ferent from the corresponding bulk material. In this re- 
spect, research on nanocatalysis is inspired by the expec- 
tation that due to such a reduced dimensionality of the 
active material new catalytic behavior can arise, which 
is not scalable from bulk-like propertiesi A salient tar- 
get of such investigations are nanometer-size clusters of 
atoms, representing the three-dimensional confinement 
situation. While the potential of such clusters for cat- 
alytic applications is just being explored, it has already 
been recognized that they often exhibit the undesirable 
propensity to coarsen into larger aggregates on the sup- 
port material. This problem could be less severe for the 
case of confinement in only one dimension, i.e. when 
striving for appealing properties of nanometer or sub- 
nanometer thin films. In fact, at the catalyst surface 
such films may even form automatically during the in- 
duction period, as a consequence of the exposure to the 
operating conditions. A prominent example of this is 
oxidation catalysis, where the oxygen-rich environment 
might oxidize the sample. Recent experimental and the- 
oretical studies noted that particularly at late transition 
metals (TMs) few or even one atomic-layer thin, so-called 
surface oxide films are formed, when these are held in a 



realistic gas-phase environment under catalytic or alike 
(T, p)-conditions. The vertical confinement makes these 
films distinct to surfaces of both bulk metals and bulk 
oxides, and their specific properties can have beneficial, 
as well as detrimental effects on the desired functionality. 

Appropriately tailoring and exploiting these function- 
alities requires an atomic-scale characterization and un- 
derstanding of the oxidation process, as well as of the 
resulting surface composition and stoichiometry (under 
the steady-state reaction conditions). Although oxide 
formation has been of technological and scientific inter- 
est since long, such an understanding is, however, only 
just emerging, and the underlying microscopic processes 
are unfortunately still largely unclear for everything that 
goes beyond sub-monolayer oxygen adsorbate phases on 
low-index single crystal surfaces. This holds in partic- 
ular for the transition from an oxygen adsorbate layer 
to the surface oxide. A main reason behind this lack of 
understanding is that atomic-scale investigations of TM 
surfaces on the verge of oxide formation still pose a signif- 
icant challenge for surface scientists: The need for rather 
high oxygen partial pressures (high compared to those 
possible in standard surface science experiments) and el- 
evated temperatures required to initiate the process, a 
low degree of order at the surface, and often complex, 
large unit-cell geometries name but a few of the obsta- 
cles encountered. Since this often exceeds the analytical 
capabilities of one experimental technique alone, multi- 
method approaches (e.g. combining scanning, diffraction 
and spectroscopic measurements) have become an impor- 
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tant tool in this field. In addition, first-principles calcu- 
lations, almost exclusively employing density-functional 
theory (DFT), have played a key role in this research 
area, in particular when suitably combined with concepts 
from thermodynamics or statistical mechanics. Although 
most of the arguments and understanding presented in 
this review are based on such calculations, I will not de- 
scribe the methods themselves here, but refer to excellent 
monographs (see for example Refs. and0) for DFT it- 
self, to Refs. 0and|^for applications to surfaces, and to 
Ref. 1^ for the hybrid statistical mechanics methods. 

Instead, this paper will focus on recent progress in the 
modeling and understanding of oxide formation of late 
TMs, using particularly the late id series from Ru to 
Ag to discuss trends. Section II will highlight aspects 
leading to the formation of surface oxides, especially on- 
surface adsorption and accommodation of oxygen in the 
sub-surface region. Key factors (also for the surface ox- 
ide stability itself) are the bulk-oxide heat of formation, 
diffusion energy barriers, elastic strain energies, and the 
nature of the surface-oxide-metal interface. Sometimes 
the nano- or sub-nanometer thickness of the formed ox- 
ide film is ruled by the thermodynamic conditions, but 
often the thickness is determined by kinetics, e.g. by 
the diffusion of metal or oxygen atoms (between inter- 
face and surface), or by the on-going oxidation reaction 
at the surface in a catalytic environment. Examples for 
both situations will be discussed in Section III, where we 
will also see that such thin oxide films can have a different 
structure to what is known from the stable bulk oxides. 
Setting these results into a "catalytic perspective" , I will 
conclude in Section IV, for example that the catalytic ac- 
tivity of Ru in high-pressure oxidation reactions is likely 
due to nanometer thin, but nevertheless bulk-like oxide 
patches, whereas for Pd and Ag the catalysis is more 
likely actuated by novel, surface-confined sub-nanometer 
thin oxidic structures. 



II. INITIAL OXIDATION OF 
TRANSITION-METAL SURFACES 

A. Formation of adlayers 

Oxide formation at metallic surfaces begins with the 
dissociation of the O2 molecule. This process is still not 
fully understood (see for example Ref. and references 
therein) , but for the present review we don't need to dis- 
cuss this in detail. As the end product of the dissociation 
atomic oxygen is chemisorbed at highly coordinated sur- 
face sites. I will mainly address the lowest energy surfaces 
of the late 4d TMs which have the hcp(OOOl) geometry 
for Ru, and the fcc(lll) geometry for Rh, Pd, and Ag. 
For these systems the oxygen adatoms occupy threefold 
hollow sites which are sketched in Fig. These two 
adsorbate positions, named hep and fee sites, differ only 
with respect to the second substrate layer: Below the 
hep hollow site there is a second-layer atom, and below 
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FIG. 1: Highly-coordinated adatom sites at a fcc(lll) or 
hep (0001) surface. The upper picture shows a top view of 
the surface, labeling the two threefold coordinated hollow 
sites. The lower row illustrates the local geometries of the 
three high-symmetry interstitial sites between the first and 
second substrate layer. Metal atoms are depicted as light, 
large spheres, and oxygen atoms as dark, small spheres. 

the fee hollow site there is none. In other words, these 
sites represent positions atoms would occupy, if the crys- 
tal were continued in a hep or a fee stacking sequence. 
Interestingly, on the four Ad materials covered in this pa- 
per, oxygen prefers initially those threefold hollow sites 
that correspond to the stacking type of the substrate, 
i.e., the hep hollow sites on Ru and the fee hollow sites 
on Rh, Pd, and Ag. With increasing number of oxygen 
adatoms repulsive lateral interactions between the adsor- 
bates lead to the formation of ordered superstructures, 
typically with a (2 x 2) surface unit-cell at low coverages. 
Then, depending on the material and whether higher cov- 
erage on-surface structures are formed at all, this may be 
followed by either (2 x 1) or (\/3 x V3)i?30° periodic over- 
layers, with oxygen atoms still occupying either hep or 
fee hollow sites. 

The quantity ruling the adsorption-site preference and 
the formation of ordered adlayers is the adatom binding 
energy: 

Eb{Q) = [So@surf. - £^surf. ^ {No/2)Eo,] ■ (1) 

Here, A^o is the total number of O atoms in the unit-cell 
at the considered coverage 9 (measured in monolayers, 
ML, i.e. the fraction compared to the number of first 
layer metal atoms). £^o@surf.i -E'surf., and Eq^ are the to- 
tal energies of the adsorbate system, of the correspond- 
ing clean surface before adsorption, and of the isolated 
oxygen molecule. A positive binding energy Ef,{Q) re- 
fiects that the dissociative adsorption of O2 is exother- 
mic. Equation |^ can be evaluated by DFT, and for big 
systems this is in fact the only approach that is currently 
feasible and predictive. One has to stress, however, that 
for oxygen at transition metal surfaces such calculations 
are still demanding. Without describing the technical 
challenges of performing such large-scale DFT studies, I 
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FIG. 2: Computed (DFT-GGA) binding energies ^^(e) for 
on-surface oxygen chemisorption into the most stable hollow 
sites of the basal surface of the late 4d TMs (Ru = hep site, 
Rh/Pd/Ag — fee site) and employing (2 x 2) surface unit-cells 
(the lines are only to guide the eye). The energies are given 
with respect to the theoretical O2 binding energy, cf. eq. Q , 
while the upper horizontal dotted line indicates the zero-level, 
if the experimental O2 binding energy were used (after Ref. 



only mention the basic problem: Though DFT is in prin- 
ciple exact, in practice it is not, because the exchange- 
correlation (xc) functional is (and probably ever will be) 
only known approximately. For many systems this is not 
a severe problem, but for the O2 molecule even present- 
day gradient-corrected (GGA) functionals suffer from a 
binding energy error of about 0.5 eV per O atomiSi2iili 
Though some error cancellation will occur in eq. (Q, this 
O2 problem calls for caution when judging on the exo- or 
endothermicity of oxygen adsorption or oxide formation. 
Binding energy differences between different adsorbate 
systems, for which the gas phase O2 problem cancels out 
at least partially, appear fortunately often much less af- 
fected by the xc approximation. Still, caution is advis- 
able, and if doubts exist, a regional xc correction, e.g., 
following the approach of Filippi et al}^, may be neces- 
sary. 

Evaluating eq. is most useful for identifying adsorp- 
tion sites, lateral interactions or the saturation coverage, 
and by such calculations the interplay between theory 
and experiment has indeed been most stimulating and 
synergetic in this field. The atomic scale characteriza- 
tion of higher coverage O adsorbate layers forms already 
an intriguing example for this. Although not really fully 
addressed and understood, it seems that the sticking co- 
efficient for O2 drops either dramatically after a certain 
threshold coverage at the surface is reached (Ru, Rh, 
Pd) or is very low from the beginning (Ag) . Particularly 
the sharp drop after reaching a coverage of about 0.5 ML 
at Ru(OOOl) and Rh(lll) leads to an apparent uptake 
saturation, when employing low gas exposures as typical 



for many early ultra-high vacuum (UHV) surface science 
studiesii2ii2*ii DFT calculations predicted, however, that 
much higher coverages should still be possibl9iS*iS*ii, 
and Fig. [3 shows such data for some selected adsor- 
bate configurations spanning the whole coverage range 
up to 1 MliiS. Although the decreasing binding energies 
with coverage indicate repulsive interactions between the 
adsorbed O atoms, they are still clearly exothermic up 
to the full ML for Ru (and also for Rh). As for the 
0/Ru case adsorption into other sites (e.g. below the 
surface) or other structures can be ruled out (see below), 
these results imply that energetically the formation of 
denser overlayers is possible, well beyond the apparent 
saturation at « 0.5 ML. And in fact, by using more 
oxidizing gases that create higher partial pressures (e.g. 
NO2), orders of magnitude higher O2 exposures (some- 
times millions of Langmuirs) or the more novel, direct 
in-situ studies at elevated pressures such overlayers (and 
configurations corresponding to even much higher O up- 
takes) are now routinely observed experimentally. One 
should stress, that on all four TM surfaces discussed in 
this paper, such experiments have opened a new door to 
study the oxidation beyond sub-monolayer overlayers in 
a controlled fashion, significantly improving the scientific 
understanding over the last years. 

Figure |21 also illustrates the mentioned 02-molecule 
problem. Using in eq. the experimental O2 total en- 
ergy, E02 , instead of the theoretical one, would shift the 
energy zero to the upper dotted horizontal line. This is 
a too extreme alternative though, because some xc error 
also exists in i?o@surf. and i^surf.- Thus, the true energy 
zero is probably somewhere in between the two dotted 
lines. Obviously, in particular for silver the inaccuracy 
in the molecular binding requires special attention, even 
when only discussing the sign of the binding energy 
What is, however, unambiguously evident from Fig. [Sjis 
that the oxygen-metal bond strength decreases steadily 
from Ru to Ag, and this is in fact known qualitatively 
since long. The reason behind it is clearly visible in the 
projected density of states (DOS), which is defined as^ 

A^o(e) = 5]|(</'c.|^,)P<5(e-e.)- (2) 

Here (pi are the Kohn-Sham orbitals, and (f>a is a prop- 
erly chosen localized function, in the present case most 
suitably the Ru Ad and the O 2p atomic orbitals within 
a sphere around given atomic sites. According to the 
Andcrson-Grimley-Newns model of chemisorption the in- 
teraction of the O 2p level with the narrow 4d band of the 
TM surface gives rise to bonding states at the lower edge 
of the c?-band and antibonding states at the upper edge 
of the d band. The corresponding peaks can be readily 
identified in Fig. |31 illustrating the projected DOS of a 
high coverage O adsorbate layer for the late Ad TM series. 
The Fermi level moves from the d-band center (for Ru it is 
already at a « 70% filling) to about 3.5 eV above the top 
of the d band (for the noble metal Ag), and thus the anti- 
bonding oxygen-metal states become progressively more 



DOS (arl). units) 



FIG. 3: Computed (DFT-GGA) projected density of states (PDOS), cf. eq. Q, of the clean basal surfaces of the late 4d 
TMs (gray shading), and of the 0(1 x 1) adlayer (black shading). Shown are the Nc,{e) contributions of the Ru 4ci and O 2p 
states. The decrease of the oxygen PDOS band width (from left to right) and the increased filling of antibonding states (the 
higher-energy oxygen PDOS region) is clearly visible. 



and more occupied 1^*22*21 As a consequence, the bond 
strength is weakened, explaining the steadily decreasing 
binding energies of Fig. |21 For the discussion below, it 
is finally worthwhile mentioning that despite these large 
variations of the binding energy with element (as well as 
with coverage), the oxygen- metal bond length does not 
change appreciably and stays always roughly at 2 A. 



B. Oxygen accommodation below the top metal 
layer 

It is often assumed that deposition of oxygen on metal 
surfaces gives initially rise to stable overlayers (either or- 
dered or disordered). Apparently this is often correct, 
but there is no reason that it holds in general. In fact, 
oxygen atoms may well go below the surface, either stay- 
ing in the near-surface region or dissolving into the bulk. 
Or, metal-oxide patches could form immediately. Exper- 
imentally, very little is unfortunately still known about 
this. For the basal surfaces of the late Ad transition met- 
als recent DFT-GGA calculations indicate that initially 
oxygen adatoms stay at the surface, and only when the 
coverage exceeds a certain value some adatoms will go 
sub-surfaceii2i They then remain preferentially close to 
the surface, i.e., they will not dissolve into the bulki22i2^ 
The interaction between these sub-surface oxygen atoms 
is attractive, and this situation may well already be called 
a surface oxide, or "stage I" of oxide formation. 

In principle, oxygen may occupy interstitial or substi- 
tutional sites in the crystal, but mostly due to the high 
cohesive energy of the late TMs the latter are found to 
be energetically much less favorable. Thus, they will 
not play a role in thermodynamic equilibrium, or close 
to it, and we may focus our discussion on interstitial 
sites. Directly below the topmost metal layer, i.e. in 
the immediate sub-surface region, there are three high- 
symmetry possibilities for such interstitial sites (cf. Fig. 
^ . Namely, there is an octahedral site (henceforth called 



octa) and a tetrahedral site (tetra-I) directly below the 
on-surface fee and hep site, respectively. And there is 
a second tetrahedral site (tetra-II) directly below a first 
layer metal atom. The main trends in the energetics of 
on- versus sub-surface oxygen are already discernible by 
comparing surface structures containing from one up to 
four O atoms in a model subset of (2 x 2) unit-cells. 
Of these up to four adatoms, consider that any one can 
be located either on- or below the surface in any of the 
high-symmetry sites displayed in Fig. ^ Despite the 
ostensible simplicity of this model subset, the number 
of possible combinations of on- and sub-surface oxygen 
atoms in various site and symmetry configurations is still 
quite large, as exemplified in Fig. ^by all configurations 
where one out of a total of four oxygen atoms is located 
below the surface. Interestingly, from all these structural 
possibilities the combination of on-surface O in fee sites 
and sub-surface O in tetra-I sites was in all four late M 
metals found to be energetically either most stable or 
very close to the most stable geometryiiSii^iSSiSiiS^ 

Occupation of sub-surface sites will commence at a cer- 
tain total coverage Oc, when a structure with a non-zero 
fraction of the O atoms below the surface is energeti- 
cally more favorable than the pure on-surface adsorp- 
tion (corresponding to a zero fraction of oxygen atoms 
below the surface). Figure |31 summarizes the energetics 
within the described model subset at the basal surface, 
and we see that for all four late Ad TMs on-surface ad- 
sorption is initially more favorable than sub-surface in- 
corporation. However, with increasing coverage this pic- 
ture changes as the aforedescribed repulsive lateral in- 
teractions between the on-surface oxygen atoms progres- 
sively diminish the preference for adsorption into these 
sites. As can be deduced from Fig. Othis threshold cov- 
erage Oc, is beyond IML for Ru and Rh, i.e. oxygen 
incorporation will only start after completion of a full 
0(1 X 1) overlayer, in nice agreement with experimental 
findings i^SiSi Close to = 1 ML the energy differences 
are very small in the case of Rh though. Hence, a fi- 
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FIG. 4: Top view of all possible on-surface/sub-surface site combinations at a total coverage = 1 ML with one oxygen atom 
per (2 X 2) surface unit-cell located below the surface, i.e. corresponding to a fraction of 25 % below the surface (see Fig. 
for the explanation of the different sites). The metal atoms are drawn as big spheres (white for the surface layer, gray for 
second layer), and the oxygen atoms are drawn as small spheres (black for on-surface, gray for sub-surface). Oxygen atoms in 
tetra-II sites below the first layer metal atoms are invisible in this plot and are schematically indicated by small white circles. 
Symmetry inequivalent occupation of the same kind of on- and sub-surface sites are denoted with (a) and (b) , respectively. 
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FIG. 5: Averaged binding energy Eb (using DFT-GGA and 
eq. (0) as a function of the total O coverage of which 
a certain percentage is located below the basal surface (on- 
surface in fee and sub-surface in tetra-I sites). The highest Et 
are always found for the pure on-surface chemisorption phase 
at = 0.25 ML (chosen as zero reference, black area), with 
each contour line (lighter gray areas) at 0.1 eV steps towards 
less stable Eb (from Ref. Uj). 



nite concentration of sub-surface oxygen might already 
be present at total coverages even slightly below 1 ML 
at elevated temperaturesiS4i2& For Pd, 8c lies around 
3/4 ML, and for Ag sub-surface incorporation starts al- 
most immediately (cf. Fig. [Sjl . 



This trend for sub-surface accommodation of oxygen 
atoms at the late Ad TMs can be understood by looking at 
the sub-surface geometries. In all relaxed structures the 
sub-surface oxygen-metal bond lengths are slightly larger 
than 2 A, i.e. they are similar to the on-surface 0-metal 
distances. Comparing these bond lengths with the avail- 
able geometrical space inside a frozen metal lattice, one 
finds them to be rather incompatible. For the unrelaxed 
tetra-I site of the late 4d metals the bond length would 
range from 1.65 A (Ru) to 1.80 A (Ag). Thus, oxygen 
incorporation necessarily must go together with a sub- 
stantial local expansion of the metal lattice. As shown 
by Todorova et ali^ the bond strength of sub-surface 
oxygen atoms may therefore roughly be divided into a 
stabilizing contribution due to the chemical binding and 
a destabilizing contribution from the lattice deformation 
cost. Going from Ru to Ag the first component should 
follow the same d-band-filling trend of the oxygen-metal 
bond strength that was discussed in Section IIA for on- 
surface oxygen. And also the latter component follows 
this trend, i.e., it scales roughly as the bulk cohesive en- 
ergy or the bulk modulus, which are largest for Ru and 
smallest for Ag as noted in Fig. |H| The two contributions 
(0-metal bond strength and metal-lattice distortion en- 
ergy) have opposite sign and apparently roughly com- 
pensate each other, leading to a rather similar energetic 
stability of sub-surface oxygen in all four elements?!^ 

Due to the lattice-deformation cost sub-surface oxy- 
gen turns out initially less favorable than on-surface ad- 
sorption, but with increasing coverage the repulsive lat- 
eral interactions for on-surface oxygen drive this pref- 
erence down until eventually penetration of some oxy- 
gen adatoms becomes favorable. As the on-surface bond 
strength decreases from Ru to Ag, while the sub-surface 
bond strength stays roughly constant, this crossover will 
occur progressively earlier, i.e. Qc is large for Ru and 
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FIG. 6: Trend of DFT-GGA computed bulk modulus, cohe- 
sive energy, and layer distance in (111) direction for the late 
4d TMs. 

small for Ag. 

If oxygen penetrates, the deformation cost also ratio- 
nalizes, why sub-surface incorporation was found to be 
more favorable compared to bulk dissolution. Although 
the latter occurs of course always to some extent just on 
entropy grounds, there will at any rate be an enrichment 
of oxygen atoms in the easier to relax interstitial sites 
in the immediate near-surface fringe. And one could ex- 
pect this deformation argument to hold even more for 
sub-surface sites close to even lower coordinated atoms, 
i.e. interstitials in the vicinity of point defects, steps or 
dislocations. The latter are in fact frequently believed 
to be the nucleation centers for surface oxide formation, 
and kinetic arguments like an easier O penetration are 
often put forward as explanation. Yet, we see that the 
thermodynamic deformation cost factor could also favor 
an initial oxygen accommodation close to such sites. 



C. Oxygen accumulation in the surface region and 
surface oxide formation 

Thermodynamically the bulk oxide will start to form at 
the oxygen coverage, 8*'^'^, at which the average oxygen 
binding energy, cf. eq. equals the oxide heat of 

formations^. The word "coverage" is used here in a wider 
sense, i.e. meaning the total amount of oxygen atoms in 
the near-surface region. Interestingly, the just discussed 
Qc for the initial incorporation of oxygen at sub-surface 
sites in the late Ad metals does not only follow O''^'^ in 
trend, but also in approximate valueiiS^. This suggests 
sub-surface oxygen as an intermediate, or a precursor, for 
the formation of the bulk oxide or other oxidic structures 
beyond on-surface adsorbate layers. Recent experimental 
studies seem to support this view, i.e., they indicate that 
thin oxidic structures, now mostly coined surface oxides, 
start to form roughly at coverages Q^,^2i^ii22i^i24 

For the "true" bulk oxide a higher oxygen concentra- 
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FIG. 7; Binding energies, calculated by DFT-GGA, for oxy- 
gen at Ru(OOOl). Goverages with B < 1 ML correspond to 
pure on-surface adsorption (hep or fee site). For O > 1 ML 
an 0(1 X 1) arrangement is always present at the surface, cf. 
Fig. 1^ while the remaining O is located in either of the three 
interstitial sites between the first and second substrate layer. 
From the six possible structures at each coverage, the three 
with on-surface O in hep (fee) sites are drawn with thicker 
(thinner) lines (from Ref. |23) . 



tion might be required, but for oxide nucleation and for 
the formation of oxidic islands such increased oxygen con- 
tent is only needed locally. In this respect it is remarkable 
that although for the on-surface adsorption at the late Ad 
TM basal surface the adatom-adatom interaction is re- 
pulsive (see the discussion of Fig. 12 above), an attractive 
interaction was found for sub-surface oxygen atoms in 
Ru(OOOl) and Rh(lll)2ai24. This is illustrated in Fig. 
{7\ which shows the course of the average binding energy 
of O at Ru(OOOl) upon increasing the amount of oxygen 
atoms between the first and second substrate layer up to 
1 ML sub-surface coverage, i.e., 2 ML total coverage. As 
discussed in the last section, for Ru(OOOl) oxygen pen- 
etration starts only after the full on-surface monolayer 
coverage is reached. Up to Q = IML the binding energy 
decreases steadily due to the repulsive lateral interac- 
tions. For even higher coverages, i.e. upon the ensuing 
filling of sub-surface sites, the binding energy curve ex- 
hibits an inflection. This implies overall attractive inter- 
actions that favor an accumulation of sub-surface oxy- 
gen in dense two-dimensional islands with a local (1x1) 
periodicity.-^ 

While for Ru(0001)22i and Rh(lll).?i the attractive in- 
teractions lead to dense oxygen islands with a local (1x1) 
ordering below the surface, cL Fig. [71 at Ag(lll) it 
already becomes more favorable to incorporate oxygen 
rather between deeper layers, when a local sub-surface 
oxygen coverage of ~ 1/4 ML is exceeded.^^ Apparently, 
each metal is only able to sustain a certain amount of 
oxygen per layer, and the differences between Ru and Ag 
are significant. Whereas Ru(OOOl) can accommodate a 
full oxygen monolayer on the surface plus a full oxygen 
monolayer below the top Ru layer, for Ag(lll) this is only 
25% of a monolayer on and 25% below the top Ag layer. 
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What seems to be a more general feature, however, is the 
0-metal-O trilayer type structure, that results from the 
accumulation of sub-surface oxygen below the topmost, 
oxygen- adsorbate covered metal layer. This accumula- 
tion will eventually lead to the formation of an ordered 
structure, which as already mentioned above, one could 
denote as a surface oxide, or "stage I" of oxide forma- 
tion. This is to distinguish it from sub-surface oxygen, 
which often conveys more the notion of a randomly incor- 
porated O species in the near-surface fringe. Although 
we expect from the deformation cost argument an en- 
richment of the latter compared to bulk dissolved O, one 
should stress that an experimental identification of either 
(bulk or sub-surface O) species will be very difficult: The 
concentration of bulk dissolved O is thermodynamically 
quite low in the late 4d TM crystal latticea^^'^^ and will 
thus only yield small signatures (which as a note aside 
does not mean that the total amount of dissolved O in the 
whole crystal must be small though). And even though 
possibly enriched to a higher concentration, the signals 
from sub-surface O will also be small, since now only a 
small volume of the crystal will contribute. In fact, sub- 
surface oxygen has hitherto been quite elusive with only 
very few experimental studies claiming to have provided 
direct evidence for it. What is more frequent is that spec- 
troscopic features that could not otherwise be explained 
have been interpreted as being due to a sub-surface oxy- 
gen species. 

Within the deformation cost picture, it appears in any 
case rather likely that (depending on the O exposure and 
kinetic barriers) a (lateral and vertical) aggregation of 
sub-surface oxygen will directly lead to ordered surface 
oxide structures. We will see below, that these struc- 
tures are indeed oxide precursors, which may often have 
a somewhat similar stoichiometry as the bulk oxides (e.g. 
Ru02 and Ag20 for the two mentioned trilayer struc- 
tures, respectively), but the geometries can differ sig- 
nificantly. Experimentally it may be similarly involved, 
though not impossible, to identify such structures, be- 
cause also an ensuing transition to other oxidic precursors 
or eventually the "true" bulk oxide (film) may proceed 
fast. 



D. Formation of the bulk oxide 

Continued accommodation of oxygen at the surface 
will eventually yield to the phase transition towards the 
"true" bulk oxide structure. For different materials this 
happens at different oxygen load, and the pathway from 
the first formed surface oxide to the final bulk oxide may 
well comprise other oxidic precursor structures. In gen- 
eral, one has to admit that our atomic-scale knowledge 
specific to this step in oxide formation is again very shal- 
low. Particularly for the more noble TMs the low sta- 
bility of the bulk oxides requires either very high oxy- 
gen pressures and/or rather low temperatures for the ex- 
perimental preparation. At the lower temperatures the 
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FIG. 8: Side view of the detached "floating trilayer" of O at 
Ru(OOOl), corresponding to the most stable geometry with 
= 3 ML. The left half of the figure shows the schematic ge- 
ometry together with the relative layer expansions compared 
to the bulk distance and the right half the calculated electron 
density. Ru = large spheres, O = small spheres, Ru atoms 
not lying in the shown plane are whitened (from Ref. l3^ . 



growth is largely affected by kinetic limitations though, 
so that at Pd(lll) at best small PdO clusters without ap- 
parent crystalline order were hitherto reported^, while 
at Ag(lll) there exists literally no atomic-scale knowl- 
edge on the oxidation process beyond the formation of an 
ordered p(4 x 4) surface oxide (see below). This situation 
is far better for Ru(OOOl) and Rh(lll), where it is at 
least known that the oxidation process ends with the for- 
mation of crystalline rutile-structured RuO2(110)^^ and 
corundum-structured Rh2O3(0001) films'^Si, respectively. 

This detailed knowledge of both initial state (the afore- 
described 0-Ru-O trilayer) and final state (the rutile 
structured RuO2(110) film) made it possible to elabo- 
rate for the first time on an atomistic pathway for the 
structural phase transition towards the bulk oxide^^: If 
the oxygen load at the Ru(OOOl) surface is increased be- 
yond the 2 ML coverage 0-Ru-O trilayer, cf. Fig. [7| 
the surprising finding of the corresponding DFT-GGA 
study was that these additional oxygen atoms would also 
go preferentially between the first and second substrate 
layer, i.e. where already one full oxygen layer exists. The 
result is displayed in Fig. |S| It shows that a significant 
relaxation takes place such that the formed 0-Ru-O tri- 
layer is essentially floating on an O /Ru adsorbate system. 
With the trilayer decoupled, increasing the oxygen load 
further then gives rise to two O-Ru-O trilayers and so 
forth. Still, this is not yet a know Ru-oxide structure, 
neither with respect to the energy, nor the geometry. 

As the trilayer on top has apparently only negligible 
influence on the underlying metal, a continued exposure 
to oxygen is likely to produce sequentially more and more 
of these trilayers, which one after the other will become 
decoupled from the metal substrate. This way a loosely 
coupled stack of trilayers will successively be formed. 
Comparing the geometry of these trilayers with the one 
of a rutile (110) layer, it can be noticed that the former 
can be transformed into the latter by a simple accordion- 
like lateral expansion of the trilayer as explained in Fig. 
1^1 With this expansion the self-contained trilayer un- 
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FIG. 9: Atomic geometries of the 0-Ru-O trilayer shown in 
Fig. |S1 (left-hand side) and of the final rutile RuO2(110) ox- 
ide structure (right-hand side). Top panel: top view, indicat- 
ing the largely different sizes of the surface unit-cells. Mid- 
dle panel: perspective view. The bulk-like rutile structure is 
achieved by simply stretching the trilayer in direction of the 
arrows, keeping the length of the drawn O-O bonds (white 
and black lines) rigid. This way an alternative sequence of 
fully coordinated Ru atoms (inside the tilted O octahedra, 
black bonds) and coordinatively unsaturated (cus) Ru atoms 
(surrounded by only four O atoms, white bonds) is created. 
Bottom panel: coupling of two trilayers vs. that of rutile 
layers. Whereas the self-contained trilayers hardly bind, new 
bonds can be formed in the RuO2(110) structure. The bridg- 
ing O atoms of the full octahedra (black bonds) form the apex 
atoms of the new octahedra (white bonds) around the former 
cus Ru atoms. Note, that the 'bonds' drawn do not refer to 
chemical bonds, but are merely used to guide the eye (from 
Ref. m. 



the effect of strain. As already pointed out before, oxygen 
incorporation will always lead to a substantial local ex- 
pansion of the crystal lattice of late Ad TMs, and known 
bulk oxide structures of the latter are accordingly much 
more open compared to the close-packed bulk metals. In 
the direction perpendicular to the surface, nothing op- 
poses such an expansion and the corresponding vertical 
relaxations upon O incorporation can be quite substan- 
tial. Lateral relaxations, on the other hand, will con- 
flict with the binding of the formed oxidic fringe to the 
underlying metal lattice. In the formation of (surface) 
oxides, the system therefore has to find an optimum en- 
ergy compromise between elastic strain relaxation and 
good coupling to the underlying substrate. This may 
result in commensurable or incommensurable overlayers, 
and both types have already been observed experimen- 
tally on the late M basal surfaces i^^'^^i??'?** In addition, 
strain may be easier relaxed at more open sites like steps, 
pointing again at the relevance of the latter for the ox- 
idation process. Apart from the strain energy, another 
influential factor is the thermodynamic driving force to 
form an oxide. The latter will roughly scale with the 
bulk oxide heat of formation, which steadily decreases 
from Ru02 to Ag20. The concomitant, lower stability of 
oxidic structures per se, together with easier to deform 
metal lattices, might then render the coupling to the un- 
derlying substrate more important towards the more no- 
ble metals, i.e. the nature of the oxide-metal interface 
plays an increasing role. In this respect it is intriguing 
to notice that particularly at Pd(lll) and Ag(lll) pscu- 
domorphic surface oxide structures have been reported, 
with a structure that seems particularly suited to maxi- 
mize the interaction with the metal lattice below. 



III. IMPLICATIONS FOR OXIDATION 
CATALYSIS 

A. The role of the gas phase 



folds into a more open geometry, in which every second 
Ru atom is undercoordinated compared to the ideal six- 
fold oxygen surrounding both offered in the trilayer and 
in bulk Ru02. It is precisely these coordinatively unsat- 
urated (cus) metal atoms that may eventually lead to a 
stabilization of the rutile structure for higher film thick- 
nesses, as they can provide extra inter-layer bonds that 
can not be formed between the saturated trilayers. While 
for only one layer thick films the more compact trilayer 
offers therefore a more favorable configuration, this pref- 
erence was already found reversed at a two layer thick- 
ness, eventually inducing the transition to rutile through 
the accordion-like expansion. 

Admittedly, this idealized pathway for one specific sur- 
face can barely be able to represent the complex atomic 
rearrangements and stoichiometry changes one can ex- 
pect more generally during this oxidation step. Still, it 
contains at least one aspect that appears quite crucial: 



A salient feature emerging from studies of the kind 
briefly sketched in the first part of this review is that the 
initial oxidation of late 4cZ TMs proceeds via few atomic- 
layer thin oxidic structures, the electronic and geometric 
structure of which may deviate significantly from their 
known bulk oxide counterparts. This makes these sur- 
face oxide films distinct to surfaces of both bulk metals 
and bulk oxides, and one is led to wonder whether they 
exhibit specific new properties that can be exploited for 
a desired functionality like catalysis. In this respect, one 
has to recognize, however, that in such applications (and 
many experiments) the actual thermodynamic variable 
is not the oxygen accommodation at the surface, but the 
partial pressures and temperature in the surrounding gas 
phase. Instead of looking at the change of the state of 
the surface as e.g. a function of increased exposure (or 
uptake) as done in the preceding sections, the relevant 
question to ask is then rather what happens to the TM 
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surface given certain environmental conditions? 

Exposed to an unlimited supply of gas phase particles 
characterized by the applied pressures and temperature, 
the surface will adapt on time scales set by the kinetic 
limitations. Already these time scales could be suffi- 
ciently long to render corresponding metastable states 
interesting for applications. In fact, the classic example 
is a slow thickening of oxide films due to limitations in the 
diffusion of oxygen atoms from the surface to the oxide- 
metal interface, or in the diffusion of metal atoms from 
the interface to the surfacei^L^ Directly at the surface 
a similar bottleneck can be the penetration of oxygen, 
which might be significantly facilitated at steps and de- 
fects and thus adds to their relevance particularly for the 
formation of surface oxides also from a kinetic point of 
view. Still, although such kinetic barriers might slow the 
oxidation process down beyond noticeable time scales, in 
a pure oxygen gas phase the true thermodynamic ground 
state will eventually be reached. With respect to oxide 
formation this is often loosely equated with a completely 
oxidized sample, and it thus came as quite a surprise 
that for the more noble 4d TMs DFT-GGA indicates that 
there are in fact (r,p)-conditions, for which a nanometer- 
thin surface oxide represents the thermodynamically sta- 
ble phase. This means that in such environments a thin 
film is automatically formed at the oxidizing surface, the 
growth of which is self-limited to the potentially partic- 
ularly interesting nanometer or sub-nanometer width. 

Reactive multi-component gas phases, as typical for 
e.g. catalysis, add another twist to this situation. In this 
case, we are dealing with an open system, in which the 
supply of reactant gases comes into contact with the solid 
surface, where a chemical reaction produces a new sub- 
stance that is then transported away. The average sur- 
face structure and composition is in this case entirely de- 
termined by the kinetics of the various underlying atom- 
istic processes like adsorption, desorption, diffusion or re- 
action. Still, for given temperature and partial pressures 
of the reactants a steady-state situation may be reached 
(or even aspired for a stable catalyst performance), but 
this steady-state must not necessarily correspond to any 
thermodynamic state. Particularly for oxidation reac- 
tions, the 0-rich environment might for example con- 
stantly oxidize the catalyst surface, while the on-going 
reaction continuously consumes oxygen in such a way, 
that oxide formation never proceeds beyond a certain 
state. Similar to the thermodynamic situation described 
above, this could also lead to a self-limited thickness of 
the oxidized surface fringe. 



B. Stability of surface oxides in an oxygen 
environment 

If we initially focus on the effect of a pure oxygen gas 
phase, it is useful to first resort to a thermodynamic de- 
scription. This allows to assess the various stages in the 
oxidation process when the surface is fully equilibrated 



with its environment, and thus provides a sound basis for 
an ensuing discussion on how kinetic limitations could af- 
fect the picture. Particularly interesting in the context 
of this review is in this respect whether the thin surface 
oxides are mere kinetic precursors to bulk oxide growth 
or may actually represent a stable phase for certain gas 
phase conditions, as well as whether we can identify a 
trend over the late Ad TM series. For a surface in con- 
tact with a gas phase characterized by temperature T and 
pressure po2J * thermodynamic description takes place 
in the grand-canonical, constant (T, p)-ensemble and one 
needs to generalize eq. ^ to evaluate the Gibbs free 
energy of adsorptiori22i22ii& 

AG(/io) = (3) 

= (GoOsurf. — Gsurf. — A^OA^O " AA''mA«m) 

Here, A is the area of the surface unit-cell, and Goosurf. 
and Gsurf . are the Gibbs free energies of the oxidized and 
clean surface, respectively. fiQ and fiM are the chemi- 
cal potentials of oxygen and metal atoms, and iVo is the 
number of oxygen atoms contained in the oxidized sur- 
face structure for which AG(/io) is computed. Finally, 
AA^M is the difference in the number of metal atoms be- 
tween the reference clean surface and the oxidized surface 
structural model. 

In the difference between the two Gibbs free energies, 
the contributions due to vibrational free energy, configu- 
rational entropy and pV-tevm cancel to some extent, and 
for the oxidized late Ad TM surfaces discussed here re- 
placing the Gibbs free energy difference by the difference 
of only the leading total energy terms appears to provide 
a reasonable approximation.'*^''*- This simplifies eq. Q 
to 

AG(AAio) ~ (4) 
« ^ [A^o£^b(0@surf.) + AoA/io -f AA^mMm] , 

where i?b(0@surf.) is the average binding energy of 
oxygen in the oxidized surface, evaluated via eq. J^l, 
and where the oxygen chemical potential is now mea- 
sured with respect to the total energy of an isolated O2 
molecule as reference zero, A/io = Mo — £'02 /2- 

Eq. 10} has a rather intuitive structure: Forming 
the oxidized surface by accomodating A^o oxygen atoms 
yields an energy gain of Afo-£'b(0@surf.), that is opposed 
by the cost of taking these O atoms out of a reservoir, 
hence A'oA/io- The equivalent term AAm/^m comes only 
into play for oxidized surfaces, where the total number of 
metal atoms is different to the one of the reference clean 
metal surface, and represents then the cost of transfer- 
ring the corresponding number of metal atoms to or from 
a reservoir with chemical potential ^m- For the metal 
atoms, this reservoir is for oxidized surfaces often con- 
veniently chosen to be the metal bulk, with which the 
surface is assumed to be in equilibriumiSi^2ii^ The oxy- 
gen atom reservoir, on the other hand, is given by the 
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FIG. 10: Computed DFT-GGA Gibbs free energy of adsorption for several ordered on-surface adsorbate structures, surface 
oxides and the bulk oxide. The dependence on A^o is also translated into pressure scales at T = 300 K and T — 600 K for 
clarity. In the bottom of the figure, the "material type" which is stable in the corresponding range of oxygen chemical potential 
is listed and indicated by the shaded regions. Left: Ru(OOOl). Compared are 0(2 x 2), 0(2 x 1), 30(2 x 2) and 0(1 x 1) 
adsorbate phases with O in hep hollow sites, cf. Fig. |21 with the two 0-Ru-O trilayer structures discussed in connection with 
Figs. [7|and|Sl Right: Pd(lll). Compared are 0(2 x 2) and 0(2 x 1) adsorbate phases with O in fee sites, cf. Fig. |5| with the 
so-called "(\/6 x \/6)" surface oxide shown in Fig. Illl fthe binding energies used to construct these graphs are taken from Refs. 



surrounding gas phase, allowing to relate A/iq to tem- 
perature T and pressure po^ via ideal gas laws42i For 
a given Afj,o{T,po2) in the environment, the thermo- 
dynamically stable phase maximizes AG(A/io), and eq. 
0) can be used to compare structural models for pos- 
sible oxidation states of the surface, e.g. oxygen adlay- 
ers with different coverage versus surface or bulk oxides. 
Since i?b(0@surf.) from eq. enters into eq. I@J, all 
the caveats about DFT calculations of the oxygen bind- 
ing energy mentioned in the beginning hold here as well. 
This uncertainty is further aggravated by the additional 
approximation due to the neglected free energy contribu- 
tions, and dictates utmost care in the interpretation of 
the obtained results. 

When explicitly comparing a set of structures in the 
sketched atomistic thermodynamics approach, surfaces 
with higher O content will become more favorable the 
more 0-rich the environment. In the limit of low oxy- 
gen pressures the clean surface will obviously be most 
stable (A/io — oo, cf. eq. Q), but this will eventu- 
ally change with increasing chemical potential in the gas 
phase. The more oxygen is accommodated in a structure, 
the steeper the slope of the corresponding AG'(A/io) 
curve with A/^o- In the limiting case of the bulk ox- 
ide with its infinite number of O atoms this will finally 
yield a vertical line at a value of the oxygen chemical po- 
tential that equals the bulk oxide heat of formation per 
O atomi^ft For any higher value of A/io, the bulk oxide 



represents then the thermodynamically stable phase. 

This overall structure is nicely visible in Fig. ^| which 
compares several of the above discussed oxidation stages 
at the Ru(OOOl) surface. While for Afio < -3.1 eV the 
clean surface is most stable, four different O adsorbate 
phases with O in hep sites and increasing coverages from 
1/4 ML up to 1 ML become progressively more stable for 
more 0-rich environments. Above A/iq < — 1.7eV, bulk 
Ru02 represents the most stable phase. The stability 
range of this bulk oxide is therewith so large that it in 
fact covers all realistically attainable pressures for tem- 
peratures below T ^ 600 K, as illustrated in Fig. ^1 
In particular, one also notices that the 0-Ru-O trilayer 
structures discussed above as possible intermediates be- 
fore the bulk oxide formation are never thermodynamic 
phases, suggesting that at best they could only be stabi- 
lized by kinetic limitations in the ensuing phase transi- 
tion to the bulk oxide. 

A qualitatively similar picture has recently been re- 
ported for Rh(lll)5i, where after adsorbate phases with 
O in fee sites, bulk Rh203 becomes most stable at higher 
oxygen chemical potentials. Due to the already decreased 
heat of formation per O atom of Rh203 compared to 
Ru02, the stability range of the bulk oxide is smaller, but 
starts in contrast to the situation at Ru(OOOl) neverthe- 
less before higher coverage adsorbate phases correspond- 
ing to ~ 0.5 ML could become favorable. This limit is 
only approximate though, since one should keep in mind 
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FIG. 11: Top- and side view of the pseudo-commensurate, 
so-called "{y/E x \/6)" surface oxide structure determined on 
Pd(lll). The atomic positions do not seem to resemble the 
bulk PdO arrangement, but exhibit roughly a O-Pd-0 trilayer 
type structure (from Ref. iSSl . 



that in such atomistic thermodynamics plots only the 
stability of a given set of structures is compared, and 
there could well be yet uncharacterized, higher coverage 
structures with a sufficiently high average binding energy 
to show up as a new stable phase. Still, bulk oxide forma- 
tion appears at this surface thermodynamically preferred 
already before a full ML 0(1x1) adsorption phase, which 
is consistent with the trend towards easier O accommo- 
dation from Ru to Ag discussed in the first part of this 
review. Although a 0-Rh-O trilayer structure equivalent 
to the one discussed for Ru(OOOl) above results similarly 
as only a kinetic precursor in DFT-GGA, it is interesting 
to note that at Rh(lll) this structure could recently be 
experimentally prepared and characterized close to step 
edges 

Proceeding to Pd(lll) one observes in most aspects a 
continuation of the identified trends. As apparent in Fig. 
1101 the again reduced thermal stability of the PdO bulk 
oxide shows up by a further decreased stability range, 
which commences nevertheless already at an O chemi- 
cal potential below the one required to stabilize higher 
on-surface adsorbate coverages with 6 ~ 1/4 ML. On the 
other hand, there is also a qualitatively new feature, since 
now the experimentally observed, so-called "(-\/6 x -v/G)" 
pseudo-commensurate surface oxide structure displayed 
in Fig. Illl is not only a kinetic precursor, but represents 
the stable phase for an intermediate range of A/io-^'^This 
is therefore similar to the reported findings at Ag(lll), 
of a p(4 X 4) surface oxide that is thermodynamically sta- 
ble over a wide range of O chemical potentials, bounded 
on the one end by the stability range of a low cover- 
age adsorbate phase and on the other end by the low- 
stability AgzO bulk oxide:SiS Both the "(v^ x a/6)" 
surface oxide on Pd(lll) and the p(4 x 4) surface oxide 
on Ag(lll) seem to have a highly complex, 0-metal-O 



trilayered atomic structure that does not resemble the 
corresponding bulk oxides, cf. Fig. II It^St^ Instead, these 
structures appear particularly suited to achieve a good 
coupling to the underlying substrate, and it is this ex- 
tra contribution that enhances their stability over a wide 
(T, p)-range beyond the one of the known bulk oxides. 

Summarizing the trend, we therefore see a changing 
role of thin surface oxide films in the oxidation process 
over the late M TM series. For Ru and Rh the present 
data suggests that such precursor structures to bulk oxide 
formation can only be stabilized by kinetic limitations. 
At the more noble metals Pd and Ag they instead repre- 
sent a thermodynamically stable phase over an increasing 
range of gas phase conditions. In corresponding oxygen 
environments, these nanometer-thin films will eventually 
form on time scales set by possible kinetic limitations, but 
never grow thicker. Due to this finite thickness, the cou- 
pling at the oxide-metal interface and an atomic structure 
that can be quite different to the one of the known bulk 
oxides, one might suspect new properties that are distinct 
to those of surfaces of both bulk metals and bulk oxides, 
and could thus be of potential interest for applications. 
For the specific functionality in oxidation catalysis, one 
needs, however, also to assess the role played by the sec- 
ond reactant. 



C. "Constrained equilibrium" 

Although the average surface structure and composi- 
tion is in catalysis even under steady-state conditions de- 
termined by the kinetics of the on-going elementary pro- 
cesses, it has been suggested that a thermodynamic ap- 
proach with certain constraints could still provide a first 
insight into the effect of the reactive environment on the 
state of the catalyst smfaceM^ A straightforward gen- 
eralization of the above described atomistic thermody- 
namics concept to a multi-component environment would 
be to account in eq. for a corresponding number of 
reservoirs described by chemical potentials A/i^ , where i 
denotes the different gas phase species. Modeling the cat- 
alyst surface in the reactive environment, each reactant is 
thus represented by a corresponding reservoir, with which 
the surface is in contact. However, in a full thermody- 
namic description all reservoirs would also be mutually in 
equilibrium. Under catalytically interesting environmen- 
tal conditions, this yields directly the products as most 
stable environment, and is thus inadequate to treat the 
effect of exposure to the reactant gas phase. 

Since it is precisely the high free energy barrier for 
the direct gas phase reaction that requires the use of 
the catalyst in the first place, it is appealing to instead 
consider the surface in contact with separate, indepen- 
dent reactant reservoirs. The surface would then be fully 
equilibrated with each reactant species, but the latter 
are not equilibrated with each other. Although neglect- 
ing the direct gas phase reaction is well justified, such 
a "constrained equilibrium" approach is still highly ap- 
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FIG. 12: Stability regions of bulk oxides of the late Ad TMs in 
(A/io, A/ico)-space, as computed with DFT-GGA. Addition- 
ally, pressure scales are drawn for T = 300 K and T = 600 K. 
The white line indicates environmental conditions correspond- 
ing to P02 = Pco ~ 1 atm and 300 K < T < 600 K, and 
environments relevant for CO oxidation catalysis over these 
metals would correspond to the near vicinity of this line. 
More specifically, the small white area indicates the {T,p)- 
conditions employed in a recent experimental study by Hen- 
driksen, Bobaru and Frenken on Pd(lOO) (s. text)»^ (from 
Ref. E^. 



proximate, as it effectively also neglects the on-going cat- 
alytic product formation at the surface itself, i.e. in other 
words exactly what one is in principle out to study. In 
this respect one should stress that the whole concept of 
a "constrained equilibrium" is not designed to yield pre- 
cise answers to the microscopic state and reactivity of 
the catalyst surface. For that, the very dynamic behav- 
ior must be modeled by statistical mechanics, as will be 
illustrated below. Yet, with atomistic thermodynamics 
the large scale behavior can be traced out first, possibly 
identifying those (T, pi)-conditions where a more refined 
treatment explicitly accounting for kinetic effects is nec- 
essary. 

How this works in practice is nicely illustrated by com- 
paring the stability range of the bulk oxides of the late 
Ad TM series within the context of CO oxidation cataly- 
sis, i.e. in a gas phase formed of oxygen and C04£ The 
progressively decreasing heats of formation of the more 
noble metal oxides already discussed in the last section 
are also clearly discernible in Fig. ^1 which shows the 
oxide stability range as a function of the two chemical 
potentials in the gas phase. At very low CO chemical 
potential, the situation in a pure O2 surrounding is recov- 



ered, and the stability range of Ru02 and PdO coincides 
with the corresponding bulk oxide shaded regions in Fig. 
mil With increasing CO content in the gas phase, this 
picture changes and the stability range gets eventually re- 
duced due to CO-induced decomposition. Although un- 
certainties in the modeling of this decomposition trans- 
late presently into rather large error bars in the respective 
diagonal limits in Fig. ^1"'" a clear trend still emerges: 
Only the stability range of bulk RUO2 extends well into 
environmental conditions representative of technological 
CO oxidation catalysis, roughly indicated in Fig. ^jas a 
white line. For Rh on the other hand the stability range 
of its oxide has already decreased significantly. Even if 
the real system was in a situation close to thermodynamic 
equilibrium, either metal or bulk oxide could prevail un- 
der ambient conditions, depending on whether the exact 
partial pressures and temperatures in the gas phase cor- 
respond to a situation above or below the stability line 
(which is at present impossible to assess due to the large 
uncertainty in the calculated position of this line). The 
bulk oxides of Pd and of Ag are finally already so unsta- 
ble that they are unlikely to play a role in CO oxidation 
catalysis. The stability of Ag20 in fact being so low, that 
it does not even show up anymore within the chemical 
potential range plotted in Fig. ^1 

It is, however, important to realize that already quite 
thin oxide films at the surface would be sufficient to in- 
duce a changed catalytic functionality. In this respect, 
just considering the stability of bulk oxides as represen- 
tatives of thick, bulk-like oxide overlayers is not enough. 
Particularly at the more noble metals, we discussed al- 
ready the extended stability range of the few atomic- 
layer thin surface oxide structures in a pure O2 gas 
phase, which could similarly prevail in a reactive multi- 
component environment. While the results of Fig. ^| in- 
dicate that thick, bulk-like PdO and Ag2 overlayers will 
hardly play a role in CO oxidation catalysis, the picture 
could be quite different for the nanometer thin surface 
oxide structures. Fig. 1131 exemplifies this with corre- 
sponding results for the case of a Pd(lOO) surface in con- 
strained equilibrium with an O2 and CO environment^^ 
Similar to the aforediscussed situation at Pd(lll), oxi- 
dation proceeds also at Pd(lOO) via a trilayer-structured 
0-Pd-O surface oxide filmS^, that was recently identified 
as essentially a strained and rumpled layer of PdO (101) 
on top of the Pd(lOO) substrataii. Afihough the (101) 
orientation is not a low-energy surface of bulk PdO,^— 
it gets stabilized in the commensurate (\/5 x \/5)R27° 
arrangement of the thin oxide film due to a strong cou- 
pling to the underlying metal substrate. In line with 
our previous discussion one would then expect an ex- 
tended stability range compared to bulk PdO, that was 
indeed subsequently found both theoretically and exper- 
imentally for a pure O2 gas phaseJS. The atomistic ther- 
modynamics results summarized in Fig. 1131 suggest fur- 
thermore that this would indeed hold equally for Pd(lOO) 
in an O2 and CO surrounding. Comparing the stability 
of a large set of on-surface (co)adsorption, surface oxide 
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FIG. 13: DFT-GGA surface phase diagram of stable surface structures at a Pd(lOO) surface in "constrained equilibrium" with 
an O2 and CO environment. The stability range of bulk PdO is the same as the one shown in Fig. 1121 Phases involving the 
{^/Ex \/Z)R27° surface oxide exhibit a largely increased stability range, that now comprises gas phase conditions representative 
of technological CO oxidation catalysis (po2 = Pco = 1 atm and 300 K < T < 600 K, marked by the black line as in Fig. I12|l 
(from Ref. 13). 



and bulk oxide structures in constrained equilibrium with 
the reactive environment^, several phases involving the 
(■\/5 X \/5)R27° surface oxide are found to be most stable 
over a wide range of (T, po2jPco)-conditions that largely 
exceed the small stability range of bulk PdO. 

The range extends in fact so much that it even just 
comprises the gas phase conditions typical for techno- 
logical CO oxidation, i.e. partial pressures of the or- 
der of 1 atm and temperatures around 300-600 K. Similar 
conditions were e.g. also chosen in a recent in-situ re- 
actor scanning tunneling microscopy (STM) experiment 
of CO oxidation at Pd(lOO) by Hendriksen, Bobaru and 
Frenken^^ cf. Figs. [Hand US With increasing 02:C0 
partial pressure ratio, the authors observed substantial 
changes in the surface morphology, which they assigned 
to the formation of a thin oxidic overlayer— . Inspecting 
Fig. ^1 it is tempting to identify this overlayer as the 
already characterized (VS x \/5)R27° surface oxide, but 
this would be rushing beyond what can be safely con- 
cluded on the basis of the presently available DFT data: 
First, surface structures that are not yet explicitly com- 
pared in the atomistic thermodynamics approach could 
still change the surface phase diagram. Second, kinetic 



effects could significantly affect the oxidation state of the 
surface. This is the case, when the on-going catalytic re- 
action consumes surface species faster than they can be 
replenished from the surrounding gas phase, thus yield- 
ing an average surface structure and composition that 
can be quite distinct to the predictions obtained within 
the constrained equilibrium concept. It could in partic- 
ular be that a fast consumption of surface oxygen atoms 
reduces the stability range of the surface oxide well be- 
low the technologically interesting environmental condi- 
tions (which are even in constrained equilibrium quite 
close to the transition to the CO-covered Pd(lOO) surface 
phases). At the moment, the only safe take-home mes- 
sage from the data presented in Fig. ^|is therefore that 
instead of the already ruled-out thick bulk- like PdO film, 
the nanometer thin (-\/5 x V5)i?27° surface oxide appears 
as a likely candidate to actuate the increased catalytic ac- 
tivity connected with the morphology changes reported 
by Hendriksen and coworkers 4^ 

Still, the proximity of the technologically relevant high- 
pressure gas phase conditions to the phase transition be- 
tween surface oxide and CO-covered Pd(lOO) in Fig. I13l is 
intriguing and stimulates further speculations. Depend- 
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ing on the exact partial pressures and temperatures, ei- 
ther the surface oxide or the Pd(lOO) metal surface could 
prevail. Even oscillations between the two phases are 
conceivable as a consequence of fluctuations, when oper- 
ating under environmental conditions very close to the 
stability limit. At finite temperatures, configurational 
entropy may furthermore give rise to a complex phase 
coexistence at the surface^^*^ which would then go in 
the direction of what was in the in-situ reactor STM ex- 
periment interpreted as a continuous consumption and 
reformation of the surface oxide even under the employed 
stationary reaction conditionsi^ Instead of the one-time 
formation of a catalytically active oxide film in the case 
of Ru, this would point at a much more dynamic state of 
the Pd surface in the reactive CO oxidation environment: 
The catalytic activity could be intimately connected to 
a reversible (local) formation and reduction of the thin 
surface oxide, possibly even giving rise to spatio-temporal 
pattern formation phenomena as often reported for CO 
oxidation at Pt-group metals4^*^2i^ In such cases, also 
the domain boundaries between surface oxide and ad- 
layers could be of particular importance for the overall 
reactivity. 



D. Kinetically limited film thickness 

The constrained equilibrium description discussed up 
to now conveys the impression that a possible oxidation 
of the catalyst surface in the 0-rich environments of ox- 
idation catalysis would rather yield bulk-like thick oxide 
films on Ru, but thin surface oxide structures on the more 
noble Ad metals. This reflects the decreasing heat of for- 
mation of the bulk oxides over the late TM series, and 
seems to suggest that it is primarily at Pd and Ag where 
oxide formation in the reactive environment could be self- 
limited to nanometer or sub-nanometer thin overlayers. 
Particularly for the case of Ru, Fig. El shows that the 
gas phase conditions typical for technological CO oxida- 
tion catalysis fall deep inside the stability regime of the 
bulk oxide, indicating that thermodynamically nothing 
should prevent a continued growth of the once formed 
oxide film. 

The formation of crystalline, bulk- like RuO2(110) dur- 
ing high-pressure CO oxidation catalysis has indeed been 
observed experimentally at Ru(OOOl), but it is interest- 
ing to note that even after long operation times the film 
thicknesses never exceeded about 20 A^L^iiS. Generally, 
one could interpret this finding as reflecting kinetic lim- 
itations to a continued growth, presumably due to slow 
diffusion of either O or Ru atoms through the formed 
film.'^''^ While this could slow the thickening of the ox- 
ide overlayer down beyond hitherto measured time scales, 
it could also be that the kinetics of the on-going catalytic 
reaction affect the surface populations in such a way, that 
no further net diffusion through the film occurs and the 
nanometer thin oxide overlayer represents in fact a ki- 
netically limited steady-state in the given reactive envi- 




FIG. 14: Top view of the RuO2(110) oxide surface explain- 
ing the location of the two prominent adsorption sites (coor- 
dinatively unsaturated, cus, and bridge, br). Also shown are 
perspective views of the four steady-state adsorption phases 
present in Fig. 1151 (Ru — light large spheres, O = dark 
medium spheres, C — white small spheres). 
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FIG. 15: Left panel: Steady state surface structures of 
RuO2(110) in an O2/CO environment obtained by first- 
principles kMC calculations at T = 600 K. In all non- white 
areas, the average site occupation is dominated (> 90%) by 
one species, and the site nomenclature indicates the species 
at the two prominent adsorption sites offered by the surface: 
br = bridge sites, cus — coordinatively unsaturated sites, cf. 
Fig. 1141 Right panel: Map of the corresponding catalytic CO 
oxidation activity measured as so-called turn-over frequencies 
(TOFs), i.e. CO2 conversion per cm^ and second: White 
areas have a TOF < 10^^cm~^s~^, and each increasing gray 
level represents one order of magnitude higher activity. The 
highest catalytic activity (black region, TOF > 10^^ cm~^s~^) 
is narrowly concentrated around the phase coexistence region, 
in which O and CO compete for both br and cus sites (from 
Ref. m. 



ronment. This would e.g. be the case, if the oxidation 
reactions at the surface consume adsorbed oxygen atoms 
at such a pace that on average O penetration into the film 
occurs not more frequently than the filling of surface sites 
with oxygen from inside the film. 

To explore this in more detail one needs to go be- 
yond the constrained equilibrium description and explic- 
itly take kinetic effects into account in the theoretical 
modeling. This can for example be achieved via kinetic 
Monte Carlo (kMC) simulations, which can essentially 
be viewed as coarse-grained molecular dynamics Simula- 
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FIG. 16: Time evolution of the site occupation by O and 
CO of the two prominent adsorption sites at the RuO2(110) 
model catalyst surface shown in Fig. 1141 The temperature 
and pressure conditions chosen (T — 600 K, pco = 20 atm, 
PO2 = 1 atm) correspond to optimum catalytic performance, 
cf. Fig. 1151 Under these conditions kinetics builds up a 
steady-state surface population in which O and CO compete 
for either site type at the surface, as reflected by the strong 
fluctuations in the site occupations (from Ref. ISal . 



tions, following the time evolution of the system on the 
basis of the rates of all individual elementary processes 
occurring in the systemA^^iS Evaluating the statistical 
interplay between this manifold of atomic-scale processes, 
kMC simulations can easily cover sufficiently long time 
scales to arrive at proper averages of quantities like sur- 
face populations, even if the system is not in equilib- 
rium, but only in a steady-state as typical for catalysis. 
Such a modeling was recently performed for the CO ox- 
idation over RuO2(110), using DFT-GGA based rates 
for 26 elementary surface processes including adsorption, 
desorption, diffusion and reaction eventsiSi Of particular 
interest for the film thickness issue is the occupation of 
surface sites with O atoms in catalytically relevant gas 
environments. In this respect, the RuO2(110) surface of- 
fers two prominent adsorption sites, commonly denoted 
as bridge (br) and coordinatively unsaturated (cus) sites 
as explained in Fig. 1141"^^ 

Figure^lshows the obtained steady-state average sur- 
face populations dX T — 600 K as a function of the gas- 
phase partial pressures, together with the simultaneously 
computed turnover frequencies (TOF, in units of formed 
CO2 per cm^ per second). The latter result from evalu- 
ating the average occurrence of the reaction events over 
long time periods as a measure of the catalytic activity, 
and are intriguingly peaked around a narrow range of 
gas phase conditions. Contrary to the situation at most 
other partial pressures, the site occupation at br and cus 
sites is then not dominated by one species, i.e. either 
O or CO (or vacant). Instead, the kinetics builds up a 
surface population in which O and CO compete for ei- 
ther site type at the surface. This competition is in fact 



nicely reflected by the huge fluctuations in the surface 
populations with time apparent in Fig. 1161 The dynam- 
ics at the surface in this "active state" is furthermore 
extremely fast, and the average time adsorbed O atoms 
stay in br or cus sites before desorbing or being reacted 
away is only of the order of fraction of milliseconds. Re- 
markably, it is particularly the concentration of O at the 
surface that turns out significantly lower than predicted 
within the constrained equilibrium approach^SiM. Even 
at the bridge sites, where O atoms bind very strongly^, 
the O occupation is only 90% and not 100% as in the nor- 
mally stable surface termination. This is quite surpris- 
ing, and indicates already how effectively the on-going 
reaction consumes surface O species. It could therefore 
well be that the surface kinetics indeed limits a contin- 
ued growth of the formed oxide film under this optimum 
catalytic performance. 



E. Surface oxidation and Sabatier principle 

The partial pressures and temperatures leading to this 
high activity state of the RuO2(110) surface in the cal- 
culations, and even the absolute TOF values under these 
conditions, agree extremely well with detailed experi- 
mental studies measuring the steady-state activity in the 
temperature range from 300-600 K and both at high- 
pressures and in UJJV'^^'^^'^^. This statement comprises 
already the emerging consensus that this Ru oxide also 
forms at the surface of Ru model catalysts in high- 
pressure environmentsSiS, i.e. the high activity that was 
originally ascribed to the Ru substrate was in fact due 
to Ru02'^^^. With this understanding a puzzle is finally 
resolved that had bothered surface scientists for quite a 
while: For a long time the high CO oxidation reactivity 
of supported Ru catalysts^?., as well as of Ru(OOOl) sin- 
gle crystals^ at high-pressures stood in sharp contrast 
to exceedingly low turnover rates observed under UHV 
conditions^, thus representing a nice example of a so- 
called "pressure gap" system^. 

Looking at the computed binding energies of oxygen at 
Ru(OOOl) shown in Fig. |21the low catalytic activity of this 
surface appears in fact quite plausible. According to the 
early Sabatier principle^^ high activity would in general 
rather be expected for intermediate bond strengths at the 
surface (stable enough to form, but weak enough to yield 
the final product). In this respect the computed medium 
binding energies at e.g. Rh(lll) and Pd(lll) in Fig. [21 
conform much better with the widespread use of these 
metals in catalytic car converters, whereas the oxygen 
bonding at Ru(OOOl) seems simply too strong. The corre- 
spondingly somehow enigmatic high catalytic activity of 
Ru(OOOl) at ambient pressures only became comprehen- 
sible when experimental studies reported the formation 
of bulk-like RuO2(110) films at the surface in such reac- 
tive environments^Si^. The subsequently established— 
and above described high activity of this oxide surface 
then not only explained the pressure gap difference (Ru 
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metal in UHV = low activity, Ru02 at high-pressures = 
high activity), but also reconciled this system with the 
Sabatier principle: The RuO2(110) surface indeed fea- 
tures oxygen species that are bound with intermediate 
bond strength comparable to the one of O at Rh(lll) 
and Pd(lll)4aia 

An equivalent, though reverse explanation has also 
been suggested for the most noble of the late Ad 
metalsiSSiSi According to Fig. |21oxygen binds to Ag(lll) 
only very weakly, which in view of the Sabatier principle 
is again difhcult to combine with the known importance 
of Ag as catalyst for several oxidation reactions. How- 
ever, the bond strength of oxygen in the p(4 x 4) surface 
oxide is substantially higher and would therefore well fall 
into the more favorable medium binding energy range. 
If the surface oxide instead of Ag metal was the actually 
active surface in the high-pressure oxidation catalysis, 
the Sabatier principle would finally be reconciled with 
the complete trend of oxygen binding over the late Ad 
series apparent in Fig. |21 Both the too strong O binding 
on Ru and the too weak O binding on Ag is alleviated 
by the formation of oxides, which exhibit a similarly in- 
termediate bond strength as oxygen at Rh or Pd. 

It thus almost appears as if the interaction of oxygen 
with TM catalysts would serve to "tune" the 0-metal 
bond strength for optimum reactivity through an appro- 
priate amount of surface oxidation. Since the thermal 
stability of bulk Ag20 is already too low and in line with 
our preceding extended stability discussion, this oxida- 
tion can at Ag at best correspond to the formation of 
a thin surface oxide though. Yet, compared to the sur- 
face oxides at Pd, even the latter is probably not stable 
enough to play an active role in CO oxidation. This can 
be different in the less reducing environments of partial 
oxidation reactions, consistent with the predominant use 
of Ag catalysts for ethylene epoxidation or formaldehyde 
synthesisSi. 

These examples underline therefore already nicely the 
potential relevance of surface oxidation in the reactive 
surroundings of oxidation catalysis. In all cases discussed 
so far, this had a beneficial effect on the overall catalytic 
performance, but detrimental effects due to the forma- 
tion of an inactive oxide film are of course equally well 
conceivable. Such a catalyst deactivation through oxi- 
dation has e.g. been frequently suspected in the experi- 
mental literature concerning Rh surfacesSSSiSi In other 
cases, the catalytic activity of oxide and metal surface 
might be comparable, so that surface oxidation would 
likely remain unnoticed unless explicitly looked for by 
experimental techniques sensitive to the atomic-scale. 



IV. CONCLUSION 

Novel experimental and theoretical approaches ad- 
dressing the effect of realistic environments on surface 
structure and composition have significantly increased 
the atomic-scale understanding of the oxide formation 




FIG. 17: Cartoon sideviews illustrating the effect of an in- 
creasingly oxygen-rich atmosphere on a late TM metal sur- 
face. Whereas the clean surface prevails in perfect vacuum 
(left), finite O2 pressures in the environment first lead to oxy- 
gen adsorption phases. Apart from some bulk dissolved oxy- 
gen, the lower deformation cost will at increasing pressures 
lead to a preferential accommodation of oxygen in the near- 
surface fringe. Thin surface oxide structures are the salient 
consequence, before eventually thickening of the oxide film 
and formation of an ordered bulk compound set in. A ther- 
modynamic or kinetic stabilization of the nanometer thin sur- 
face oxide could lead to novel functionalities, different from 
both bulk metal or bulk oxide surfaces. 



process of late transition metals over the last years. Some 
identified key features have been reviewed for the late Ad 
sequence from Ru to Ag, and are summarized schemat- 
ically again in Fig. 1171 A progressively more 0-rich gas 
phase leads first to on-surface adsorption, followed by O 
penetration into the surface at increasing O accommoda- 
tion. The incorporation takes place after the formation of 
dense oxygen overlayers in the case of Ru, but already at 
lower on-surface coverages towards the more noble met- 
als. This results from the opposing trends in decreasing 
0-metal bond strength and decreasing deformation cost 
of the metal lattice, the latter being a consequence of 
the substantial local expansion necessarily connected to 
O incorporation into the close-packed TM lattices. The 
deformation cost is also responsible for an enrichment of 
oxygen in the easier to relax immediate surface fringe, 
eventually inducing the formation of thin oxidic struc- 
tures at the surface. Although some bulk dissolution will 
always occur on entropy grounds, there is thus primarily 
a propensity to form few atomic-layer thin surface oxides, 
which mostly goes hand in hand with substantial struc- 
tural changes. The resulting geometry is then not only 
characterized by an optimum binding within the more 
open compound structure, but can also be largely af- 
fected by the coupling to the underlying metal substrate 
(i.e. the nature of the oxide-metal interface) and residual 
strain. As consequence, complex atomic arrangements 
not resembling surfaces of bulk oxides or bulk metals 
have been reported for this step in the oxidation pro- 
cess. Only after continued oxygen accommodation, the 
thickening oxide film will eventually switch to a bulk ox- 
ide structure, and the continued film growth will then be 
more and more limited by the slow diffusion of O and/or 
metal atoms between surface and interface. 

A most intriguing aspect of this surface oxidation is 
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that it proceeds via the few atomic-layer thin surface 
oxide structures. This thickness can be ruled by ki- 
netics, but at Pd and Ag it even seems that there are 
thermodynamic conditions for which the surface oxides 
represent the most stable phase. In corresponding envi- 
ronments such films are formed automatically, but the 
growth is self-limited to the potentially particularly in- 
teresting nanometer or sub-nanometer width. If stabi- 
lized in applications, the often quite distinct geometries 
of these films already suggest that new properties and 
functionalities might emerge that are not scalable from 
the known bulk materials. At present, very little to noth- 
ing is known about these properties, and one should re- 
call that not even the corresponding bulk oxides of the 
late Ad TMs have hitherto received a lot of attention in 
atomic-scale studies. This opens up a wide new field for 
investigations of this new "material class" and its func- 
tionalities in various applications. 

With respect to oxidation catalysis, oxide films may 
even form as a natural consequence of the reactive en- 
vironment, and since maybe hitherto unnoticed, they 
might then actuate some of the catalytic behavior that 
has traditionally been ascribed to the metal substrates. 
In this respect, the decreasing heat of formation leads to a 
clear trend for the bulk oxides: A bulk- like Ru02 film can 
even sustain the reducing environments of high-pressure 
total oxidation reactions, whereas the feeble Ag20 will 
hardly prevail. This is, however, not enough to exclude 
the relevance of oxidized structures for oxidation cataly- 
sis over Pd and Ag catalysts. Also due to the coupling to 
the underlying substrate, the thin surface oxides exhibit 
a much extended stability range, and it seems at the mo- 
ment at least feasible that for example CO oxidation over 



Pd is actuated by them, or their constant consumption 
and reformation. The present data suggest a similar role 
for surface-confined oxidic structures at Ag, but only in 
less reducing environments e.g. typical for partial oxi- 
dation. This would then also offer an intuitive explana- 
tion for the often reported high catalytic activity of noble 
metals for such reactions. Together with the formation 
of highly active, nanometer thin Ru02 films at the sur- 
face of Ru model catalysts, this underlines the potential 
importance of surface oxidation in the reactive environ- 
ments of oxidation catalysis over late TMs. Although 
often not yet fully conclusive, already the reviewed new 
insights into the oxidation behavior help sometimes to 
reconcile the data with more general concepts like the 
Sabatier principle. Much more often, however, they lead 
to the falldown of "established" rules and preconceptions, 
which is a prospect that makes working in this field ex- 
citing and all the worthwhile. 
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